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Abstract 


The power spectral density of the output of wind turbines provides information on the character of fluctuations in turbine output. Here both 
1-second and 1-hour samples are used to estimate the power spectrum of several wind farms. The measured output power is found to follow a 
Kolmogorov spectrum over more than four orders of magnitude, from 30s to 2.6 days. This result is in sharp contrast to the only previous study 
covering long time periods, published 50 years ago. The spectrum defines the character of fill-in power that must be provided to compensate 
for wind’s fluctuations when wind is deployed at large scale. Installing enough linear ramp rate generation (such as a gas generator) to fill in 
fast fluctuations with amplitudes of 1% of the maximum fluctuation would oversize the fill-in generation capacity by a factor of two for slower 
fluctuations, greatly increasing capital costs. A wind system that incorporates batteries, fuel cells, supercapacitors, or other fast-ramp-rate energy 
storage systems would match fluctuations much better, and can provide an economic route for deployment of energy storage systems when 


renewable portfolio standards require large amounts of intermittent renewable generating sources. 
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1. Introduction 


Wind power produced by turbines varies with time. Data with 
1-s time resolution for one 1.5 MW nameplate capacity turbine 
are shown in Fig. 1. The variability is not dramatically reduced 
when hourly samples of the output of several turbines in a wind 
farm are summed, as in Fig. 2, nor when hourly data from four 
wind farms are summed, as in Fig. 3 (the median hourly change 
is reduced by 25%). 

The output of such turbines is not random. The character of 
the variations can be examined in several ways. One method [1] 
is to construct a histogram of the step size in output over time. 
Here I extend an approach that has seen use over only a limited 
frequency range [2]: constructing the power spectrum of wind. 
Data here cover four decades of frequency, using both 1-s and 1- 
h time resolution data. This frequency range is 100 times larger 
than that of previous studies. 
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2. Data 


Hourly sampled real power output data with 1 kW power 
resolution for four wind farms were obtained for the period from 
2001 to 2004. The operator of the wind farms has requested that 
they not be identified by location. Each turbine—generator in the 
four locations was identical. The minimum distance between 
the generation facilities was 30 km; the maximum was 400 km 
(Table 1). These data were provided for the output of each wind 
farm, not for individual turbines. For certain periods, the data 
were supplied with 1 MW power resolution. All four locations 
had 1 kW resolution hourly data available for the period January 
1—June 30, 2004. Wind farms A and B had 1 kW hourly data for 
July 1, 2003—June 30, 2004. 

One-second sampled real power output data with 0.1 watt 
power resolution for each of the six turbines at wind farm A and 
each of the 10 turbines at wind farm B were obtained for a 10-day 
period in 2005. Wind farm A data for the final day ceased after 
approximately the first 21 h, leaving a total of 841,600 continu- 
ous samples (Fig. 1). Wind Farm B data were good throughout, 
giving a total of 864,000 continuous samples. The individual 
turbine outputs were then summed to give the output of each 
wind farm (Figs. 4 and 5). 
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Fig. 1. Real power output (kW) sampled with 1-s resolution for one 1.5 MW 
turbine at one wind farm (farm A) for 10 days. Negative values are due to turbine 
electrical loads. 
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Fig. 2. Real power output as a percent of nameplate capacity sampled with 1-h 
resolution for the sum of the six turbines at wind farm A for the period January 
1-June 30, 2004. 
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Fig. 3. Real power output as a percent of nameplate capacity sampled with 1-h 
resolution for the sum of the 104 turbines at wind farms A, B, C, and D from 
January 1 to June 30, 2004. 


Table 1 
Distance (km) between wind farms 

A B C 
D 320 350 400 
A 30 100 
B 90 
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Fig. 4. Real power output (kW) at 1-s resolution as the sum of the six turbines 
at wind farm A over a 10-day period. 
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Fig. 5. Real power output (kW) at 1-s resolution as the sum of the 10 turbines 
at wind farm B over a 10-day period. 


3. Geographic correlation 


The 4368 hourly data points from January to June 2004 were 
used to examine the correlation of the power output between 
each pair of wind farms. Pearson’s correlation coefficient, 


Xei — HO — ¥) 
y DiGi — ay Dili — 9)? 


was used to determine the linear correlation between the real 
power outputs of the facility pairs. More detailed techniques are 
available for other correlation analysis [3]. 

The output power of wind farms that are in relatively close 
proximity are strongly positively correlated (Table 2). 


(1) 
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Table 2 
Output power correlation between wind farms, January 1—June 30, 2004 

A B C 
D 0.46 0.35 0.36 
A 0.69 0.71 
B 0.77 
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Fig. 6. Real power output as a fraction of nameplate capacity of wind farms B 
and C for the first 4368 h of 2004. Each point represents one hourly datum. 


The highest observed correlation was between wind farms B 
and C. Each point in Fig. 6 is the real power output from 1 h in 
the data set for these two facilities. 

It is sometimes asserted that combining outputs from geo- 
graphically separated wind turbines dampens most power 
fluctuations. However, even at a distance of 400 km, there is 
a fair degree of positive correlation between wind farm D and 
the other three facilities. 


4. Spectral analysis processing 


To estimate the power spectrum (sometimes termed the power 
spectral density or PSD) of the real power output of a wind 
farm, the discrete Fourier transform of the time series of out- 
put measurements is computed. An N-point-long equal interval 
time sample of the output c(t) is used to construct the value at 
frequency domain point k, Cx [4] 


N-1 
Ce= Sie", &=0,1,...,N—1 (2) 
j=0 


The periodogram estimate of the power spectrum at fre- 
quency domain point k is then [5] 


1 2 
Po = pal 
Ph = gallel + lcw] k=1,2,.. yl (3) 


P fax = 39 


where the relationship between the frequency and frequency 
domain point k is given by 

N 
2 
and fmax is (in concert with the Nyquist sampling theorem) 
0.5 Hz for the 1-s data, and 1.4 x 1074 Hz for the 1-h data. 

The periodogram is normalized so that the sum of all points P 
is the mean squared amplitude of the time series c(t). To enable 
better comparison with other attributes of the electric power 
system, the graphs that follow display the square root of the 
periodogram, so that the units are in kW /Hz~!. 

One of the attributes of power spectrum estimation through 
periodograms is that increasing the number of time samples (N) 
does not decrease the standard deviation of the periodogram at 
any given frequency fg. In order to take advantage of a large num- 
ber of data points in a data set to reduce the variance at fg, the data 
set may be partitioned into several time segments. The Fourier 
transform of each segment is then taken and a periodogram esti- 
mate constructed. The periodograms are then averaged at each 
frequency, reducing the variance of the final estimate by the 
number of segments (and reducing the standard deviation by the 
reciprocal of the square root of the number of segments). 

In this work, eight segments are used. This has no effect on 
fmax, but increases the lowest non-zero frequency (f1) by a factor 
equal to the number of segments. Thus, the 366 days of hourly 
data treated in this manner have a minimum frequency compo- 
nent of 2.5 x 1077 Hz and a maximum frequency corresponding 
to 2h (1.4 x 1074 Hz). When this technique is used, the 10-day 
time series acquired at 1-s resolution has frequency components 
from a minimum of 9.2 x 1076 Hz to a maximum of 0.5 Hz. 

Data windowing with both a Bartlett and Welch window fol- 
lowed by segment overlapping by one-half the segment length 
[6] was tried on the 1-s resolution data, with no noticeable 
improvement in the variance. The power spectra presented below 
were estimated without windowing (i.e. using a square window) 
and without overlapping the segments. In order to ensure that the 
processing introduced no artifacts in the power spectra, the wind 
farm A 1 s data were replaced by pseudorandom data (generated 
using a Marsaglia—Zaman subtract-with-borrow algorithm [7]) 
with the same maximum and minimum as the wind farm A data. 
The resulting PSD was flat on a log—log plot, as expected for 
random data. 


k 
fk = 2 fmax > k=0,1,... (4) 


5. Power spectra 


Using the method discussed in Section 4, the power spectra 
for the wind farms A and B were estimated, using hourly sampled 
data from July 1, 2003 to June 30, 2004 and using 1-s sample 
data from a 10-day period. 

The spectrum of the output of the sum of the six turbines at 
wind farm A using the hourly sampled data (Fig. 7) and the 1-s 
interval data (Fig. 8) were combined to yield a spectrum from 
2s to 6 weeks (Fig. 9). 

The spectrum of the output of the sum of the 10 wind farm B 
turbines using the hourly sampled data and the 1-s interval data 
are combined in Fig. 10. 
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Fig. 7. Power spectrum of data sampled at 1-h resolution of the sum of the six 
turbines at wind farm A from July 1, 2003 to June 30, 2004 using 8-segment 
averaging. 
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Fig. 8. Power spectrum of data sampled at 1-s resolution of the sum of the six 


turbines at wind farm A using 8-segment averaging. 


The combined spectra of Figs. 9 and 10 are characterized 
by four different regions. At frequencies between 2 x 1076 
and 4 x 107? Hz, the double logarithm plot of the spectrum 
is linear (i.e. exponential in frequency). At frequencies above 
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Fig. 9. Combined wind farm A power spectra from Figs. 7 and 8. 
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Fig. 10. Combined power spectrum of data sampled at 1-h and 1-s resolution of 
the sum of the 10 turbines at wind farm B using 8-segment averaging. 
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Fig. 11. The combined wind farm A spectra of Fig. 9 (grey points) and a f-*/ 
spectrum (solid line). 


~5 x 107? Hz, the physical and electrical inertia of the turbines 
appear to act as low-pass filters. Frequencies between 107! and 
5 x 107! Hz show the noise floor of the output power sensors. 
Peaks in the spectrum in this region are generally due to the tur- 
bine blade passing frequencies [8]. At frequencies below roughly 
2 x 107 Hz the maximum output of the turbines provides an 
upper limit to the spectrum. 

The linear region of the power spectrum plot is well fit by an 
exponential function in frequency of the form f~’? (Fig. 11). 

The combined output of all four wind farms also follows the 
same form, over the frequency range sampled by the hourly data 
(Fig. 12). 


6. Kolmogorov spectrum 


There are good theoretical reasons to believe that the power 
spectrum follows the form found in Section 5. The Reynolds 
number (ratio of inertial to viscous force) of the Earth’s atmo- 
sphere can be ~10°. For sufficiently large Reynolds number, 
viscous effects can be neglected, and the energy dissipation in 
wind is dependent on only inertial forces. In this regime, for 
an incompressible fluid, Kolmogorov’s second hypothesis of 
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Fig. 12. The hourly spectra of the combined output of all four wind farms shown 
in Fig. 3 (grey points) and a f->? spectrum (solid line). Note the scale change 
on the ordinate from previous figures. 


similarity [9] predicts that the power spectrum should vary as 
f-5/3. 

The first experimental verification of Kolmogorov’s spectrum 
was in an ocean tidal channel near Vancouver, over two decades 
of frequency [10]. Similar spectra have been observed in the 
atmosphere [11]. 

Previous published power spectra of wind generator power 
[2] have linear regions of the power spectrum plot covering one 
or two decades of frequency. These studies have not provided a 
comparison between their data and the Kolmogorov spectrum. 

Analysis of the wind power spectrum has been hampered by a 
power spectrum of wind speed from 1.9 x 1077 to 0.25 Hz pub- 
lished in 1957 [12] that has been reprinted in a recent handbook 
[13] and review paper [14]. This spectrum has a pronounced 
“spectral gap” between about 3 x 1075 and 7 x 107°? Hz with 
very little energy. The data analyzed here show no such gap, and 
an area-preserving plot (frequency times PSD against the log of 
frequency) shows a smooth behavior through this region. This 
disagreement may arise from the way in that the older spectra 
were measured (at several different altitudes) or in the way those 
data from several time periods were combined. 

The measured output power in the present study follows a 
Kolmogorov spectrum over more than four orders of magni- 
tude in frequency, from 30 s to 2.6 days. The actual atmospheric 
behavior may extend further, since these data are constrained 
by the turbine inertia at high frequencies and by the maximum 
output of the turbines at low frequencies. 


7. Discussion 


As shown in Sections 5 and 6, the output power of individual 
turbines or wind farms can effectively be modeled by a f—> 
relation for a broad range of frequencies, extending from the 
turbine inertial limit at high frequency to the wind farm power 
output limit at low frequency. 

The only previous use of the Kolmogorov spectrum for mod- 
eling of wind power systems has been for flicker analysis at 
frequencies above 10-2 Hz [8]. However, caution should be used 
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Fig. 13. The spectrum of load in one control area near the turbines used for 
this study for 180 days at 14-s time resolution using 8-segment averaging, and 
a Kolmogorov spectrum (f~~?) displaced upwards from the best fit amplitude 
for clarity. 


at frequencies grater than roughly 5 x 107? Hz for turbines such 
as those measured here, because their inertia has the effect of 
applying a low pass filter to the wind spectrum (Figs. 9-12). 

Wind at small scale is sometimes treated by grid operators as 
negative load. To examine the validity of this empirical practice, 
load data from a control area near the turbines used in this study 
were obtained with 14 s resolution for the first 180 days of 2004. 
The PSD (Fig. 13) estimated in the same way as for the wind 
power data above shows the expected peaks at 1 day and its 
harmonics, as well as the weekly peak. There is aregion, between 
1h and 2.5 min, in which the spectrum of this region’s load is fit 
well by a Kolmogorov spectrum. In this interval, the practice of 
treating wind power fluctuations as negative load appears to be 
justified. Note that this is not the same thing as saying that load 
fluctuations cancel wind fluctuations, as is sometimes stated, 
since the two would have to be both of the same magnitude and 
anti-correlated for the assertion to be valid. 


8. Conclusion 


Since wind is an intermittent resource, it must be matched 
with fill-in power sources from storage or generation if the power 
output of wind farms on a grid are correlated. 

If wind at very large scale were to be used, what ramp rate 
characteristics would be required by the fill-in energy? The 
ensemble of generators, energy storage, and demand response 
used should have a power spectrum that matches that of the wind 
generators, the f~% Kolmogorov spectrum. 

A line connecting the 24h load peak and its harmonics 
(Fig. 13) falls off as f2. This rapid decrease in power at high 
frequency allows the use of slow-ramping generators to match 
the fluctuations. However, the f~ curve for wind power means 
that high frequency fluctuations contain considerable power. 

A linear ramp rate generator is not the optimum match for 
wind. The power spectral density for such a generator has the 
form f~*. Suppose that such a generator were to be sized to 
compensate for the variations in wind power at 1% of the max- 
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imum variations observed. This point is reached at roughly 
2.5 x 1074 Hz for the systems shown in Figs. 9 and 10. Then 
(because the slope of the generator’s PSD is steeper than the 
Kolmogorov spectrum) the generator would be sized to com- 
pensate for fluctuations about twice as large as the maximum 
actually observed at low frequencies. 

A more efficient match between wind’s fluctuations and fill-in 
power can be made by noting that the source is required to match 
fluctuations at high frequency, but at greatly reduced amplitude 
compared to the match required at low frequency. 

Rather than size a linear ramp rate generator (often a natural 
gas generator in current practice) to match fluctuations in the 
output power of wind farms, a more efficient solution is to match 
wind with an ensemble of generators and energy storage to match 
the Kolmogorov spectrum of wind power. Fast devices including 
batteries, fuel cells, or supercapacitors, with relatively low power 
would match the short-period fluctuations, while slower ramp 
rate sources would match the longer period, higher amplitude 
fluctuations. The capital savings from not over-building linear 
ramp rate matching generation by a factor of two can be large 
for wind at large scale. 

The power spectra presented here can be used to evaluate 
the suitability of battery, fuel cell, and supercapacitor systems 
for use in large-scale wind systems. The advent of renew- 
ables portfolio standards in many of the United States provides 
an opportunity to deploy well-matched battery, fuel cell, and 
supercapacitor systems with wind to provide a technical and 
economically efficient solution to wind’s fluctuation. 


Acknowledgements 
This work was supported in part by the Alfred P. Sloan Foun- 


dation and the Electric Power Research Institute under grants 
to the Carnegie Mellon Electricity Industry Center, and in part 


by the U.S. National Science Foundation under awards CCR- 
0325892, SES-0345798, and CNS-0428404. The author thanks 
Joseph DeCarolis, M. Granger Morgan and José Moura for help- 
ful discussions. 


References 


[1] See for example, Y. Wan, Wind Power Plant Behaviors: Analy- 
ses Of Long-Term Wind Power Data, National Renewables Energy 
Laboratory Technical Report NREL/TP-500-36551, 2004. Available: 
http://www.nrel.gov/docs/fy04osti/3655 1 .pdf; 

Y. Wan, D. Bucaneg, J. Solar Energy Eng. 124 (4) (2002) 427-431. 

[2] P. Sørensen, A.D. Hansen, P. Andre, C. Rosas, J. Wind Eng. Ind. Aerodyn. 
90 (12-15) (2002) 1381-1402; 

G. McNerney, R. Richardson, IEEE Trans. Energy Conv. 7 (4) (1992) 
644-647; 
T. Thiringer, IEEE Trans. Energy Conv. 11 (3) (1996) 601-606. 

[3] P. Sørensen, A.D. Hansen, P. Andre, C. Rosas, J. Wind Eng. Ind. Aerodyn. 
90 (12-15) (2002) 1381-1402. 

[4] R. Bracewell, The Fourier Transform and its Applications, McGraw-Hill, 
New York, 1965. 

[5] W.H. Press, S.A. Teukolsky, W.T. Vetterling, B.P. Flannery, Numerical 
Recipes in FORTRAN: The Art of Scientific Computing, second ed., Cam- 
bridge University Press, Cambridge, 1992, p. 544. 

[6] W.H. Press, S.A. Teukolsky, W.T. Vetterling, B.P. Flannery, Numerical 
Recipes in FORTRAN: The Art of Scientific Computing, second ed., Cam- 
bridge University Press, Cambridge, 1992, p. 549. 

[7] G. Marsaglia, A. Zaman, Ann. Appl. Probab. 1 (3) (1991) 462-480. 

[8] C.V. Moreno, H.A. Duarte, J.U. Garcia, IEEE Trans. Energy Conv. 17 (2) 
(2002) 267-272. 

[9] A.N. Kolmogorov, Dokl. Akad. Nauk. SSSR 30 (1941) 301-305; 
Reprinted R. Soc. Lond. Proc.: Math. Phys. Sci. 434 (1890) (1991) 9-13. 

[10] H.L. Grant, R.W. Stewart, A. Moilliet, J. Fluid Mech. 12 (1961) 241-263. 
[11] A.S. Monin, Phys. Fluids 10 (9P2) (1967) S31-S37; 
J.C. Kaimal, Boundary Layer Meteorol. 4 (1973) 289-309. 
[12] I. Van der Hoven, J. Meteorol. 14 (1957) 160-164. 
[13] T. Burton, D. Sharpe, N. Jenkins, E. Bossanyi, Wind Energy Handbook, 
John Wiley & Sons, Chichester, 2001 (Chapter 2). 
[14] E.J. Fordham, Wind Eng. 9 (1985) 95-133. 


